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We show that all central rapidity hadron yields measured in heavy ion collisions at LHC are 
well described by the chemical non-equilibrium statistical hadronization model (SHM), where the 
chemically equilibrated QGP source breaks up directly into hadrons. SHM parameters are obtained 
as a function of geometric centrality of colliding ions, and we compare LHC with RHIC62 results. We 
predict yields of unobserved hadrons and address anti-matter production. The physical properties 
of the quark-gluon plasma fireball particle source show universality of hadronization conditions at 
LHC and RHIC. 
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I. INTRODUCTION AND MOTIVATION 



Our interest in the multi-particle production process 
in ultra-rclativistic heavy ion collisions originates in the 
understanding that the transverse momentum integrated 
rapidity distributions are insensitive to the very diffi- 
cult to fully characterize transverse evolution dynamics 
of the hot fireball source [l|. A successful description 
of central rapidity particle yields in a single freeze-out 
model 0, Q will be used here to characterize the proper- 
ties of the hadronizing quark-gluon plasma (QGP) fire- 
ball. The QGP breakup, as modeled within the statis- 
tical hadronization model (SHM), assumes equal reac- 
tion strength in all hadron particle production channels. 
Therefore, the phase space volume determines the hadron 
yields. SHM has been described extensively before and we 
refer the reader to SHARE manuals Q for both, further 
theoretical details, and numerical methods. Here, we ap- 
ply SHM to study particle production in Pb-Pb collisions 
at y/SNN — 2.76 TeV (LHC2760), a new energy domain 
an order of magnitude higher than previously explored 
at RHIC200. 

We begin by demonstrating that the chemical non- 
equilibrium SHM variant we have developed in the past 
15 years [1-0| resolves the myth that SHM does not work 
at LHC. We find that chemical non-equilibrium SHM de- 
scribes the experimental LHC-ion data with high accu- 
racy. In the chemical non-equilibrium SHM approach, 
we allow quark pair yield parameter 7g also for light 
quarks. We validate our numerical approach by show- 
ing correspondence of chemical equilibrium SHM results 
with other fits that announced failure of SHM d, Q. 
This demonstrates that several SHM programs, which 
had years time to mature and evolve, are compatible 
in their data tables of hadronic resonance mass spectra 
and decay patterns. However, only our extended SHARE- 
code includes advanced features, such as chemical non- 
equilibrium of all quark yields, differentiation of up and 



down quarks, evaluation of fireball physical properties, 
and the capability to constrain the fit by imposing phys- 
ical properties on the particle source. 

To demonstrate that our chemical non-equilibrium 
SHM works at LHC, we show in the left panel of figure [1] 
our fit to the 0-20% centrality data, shown in the second 
column of table HI recently presented and studied by the 
experiment ALICE Si, 9]. Only in this one instance, we 
consider the relatively wide geometric trigger of 0-20% 
centrality so as to compare directly to the earlier anal- 
ysis effort. As can be seen in this left panel of figure [TJ 
our non-equilibrium SHM approach describes these data 
with x^/ndf = 9.5/9 ~ 1. We see, in top figure [T] in- 
sert, that the chemical equilibrium SHM works poorly, 
X^/ndi = 64/11 ~ 6. This simple SHM disagrees at LHC 
across many particle yields but the greatest issue is with 
the p/tt — 0.046± 0.003 ratio 9], a point we will study in 
more detail in subsection IHI Bl This failure of the infa- 
mous chemical equilibrium SHM is the origin of the myth 
that SHM does not work at LHC 0, The term SHM 
is much wider than the simplified chemical equilibrium 
SHM variant. 

The chemical non-equilibrium results for the 5-10% 
centrality bin (containing interpolated data, open sym- 
bols) is shown for comparison in the right panel of fig- 
ure [H The fit has the same set of particles as the 0-20% 
centrality bin, however, we must fit here three ratios for 
which data are directly (or by interpolation) available, 
and we use a more recent set of proton, pion and kaon 
data. Definition of the model and some technical details 
about how we obtain results seen figure [T] follow below; 
the fitted data are shown in the fourth column of table H] 
The top section of the right panel of figure [T] shows the 
SHM parameters and for all three variants. Compar- 
ing the SHM parameters on left and right-hand of figure[T] 
we see a large change in V expected for different centrali- 
ties. We see that use of finer centrality binning and more 
mature data sample reduces for all SHM variants. 
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FIG. 1: (color online) The non-equilibrium SHM fit to LHC-ALICE experimental data; left panel shows 0-20% centrality, 
right panel shows 1/4 of this range, 5-10%. The upper parts show the ratio of SHM to experimental data for non-equilibrium, 
semi-equilibrium and equilibrium variants. In that part, anti-particles are omitted for readability. As seen in the in-figure 
legends, the input data sets differ for the two cases, but the set of particle types is the same. Note, in the right figure panel, 
that interpolated experimental values we use are shown with open symbols (see appendix |X] for details). 



As figure [T] shows and we discuss below in detail, the 
chemical non-equilibrium SHM works perfectly at LHC, 
resulting in a high confidence level. This could be pre- 
dicted considering prior SPS and RHIC data analysis [10|- 
[T2I , which strongly favors chemical non-equilibrium vari- 
ant of SHM. Moreover, the chemical non-equilibrium SHM 
has a dynamical physical foundation in sudden breakup 
of a QGP fireball, we are not aware of a dynamical origin 
for the simple chemical equilibrium SHM since no dy- 
namical computation of relativistic heavy ion scattering 
achieves the chemical equilibrium condition without in- 
troduction of unknown particle, cross sections, etc.. Fur- 
thermore, as we will discuss in subsection IIII C[ we obtain 
hadronization universality across a wide collision energy 
range: comparing RHIC62 with LHC2760 we show that 
the fireball source of particles is nearly identical, and con- 
sistent with chemically equilibrated QGP fireball. Given 
this result, chemical non-equilibrium SHM variant is val- 
idated across a wide energy range, while the chemical 
equilibrium SHM [l3l - [ll | is invalidated by the LHC data 
and this conclusion can be extended across different re- 



action energies as there is no reason why a model should 
work only sporadically. 

We have now shown that the chemical non-equilibrium 
is the necessary ingredient in the SHM approach to the 
process of hadronization of a QGP fireball. The non- 
equilibrium SHM was proposed when first strange hadron 
multiplicity results were interpreted more than 20 years 
ago [i^. The yield of strange hadrons indicated that 
the number of quark pairs present had to be modified by 
a factor 7^, the source of strangeness is not populating 
the final state hadrons with the yields expected from the 
hadronic chemical equilibrium, a point of view widely ac- 
cepted today. At SPS energies, for which this model was 
originally conceived, production of strangeness did not 
yet saturate the QGP phase space, that is strangeness was 
out of chemical equilibrium both in QGP fireball source 
with 7^ < 1, and thus also in the final hadronic state 
with also 7^ < 1. The distinction of QGP as initial and 
the final hadron phase space domain for 75 was also mod- 
eled poj . It is important to always remember that hadron 
phase space non-equilibrium can arise from a QGP fire- 
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ball with strangeness in chemical equilibrium, since, in 
general, the QGP and hadron phase space strangeness 
density are greatly different. Moreover, it is quite pos- 
sible that a sub-chemical equilibrium QGP, which is the 
higher density phase, produces an equilibrated hadron 
yield. This can, however, happen only accidentally and 
variation of reaction energy or collision centrality shows 
this. 

Another non-equilibrium parameter 7c, similar to 7s, 
was introduced very soon after 7^ to control the charm 
final state phase space [l^l, and it has been widely 
adopted in consideration of a strong charm yield over- 
abundance above chemical hadron gas equilibrium. Note 
that both strangeness and charm flavors are therefore 
assumed to have been produced in a separate and inde- 
pendent process before hadronization — and note fur- 
ther that each of the production mechanisms, in this 
case, is different with charm originating in first parton 
collisions and strangeness being produced in secondary 
thermalized gluon fusion reactions. At the end of QGP 
expansion, these available and independently established 
strangeness and charm particle supplies are distributed 
into available final state phase space cells, that is the 
meaning of SHM in a nutshell. 

The full chemical non-equilibrium is introduced by 
means of the parameter 7^ ^ 1. This situation arises 
when the source of hadrons disintegrates faster than the 
time necessary to re-equilibrate the yield of light quarks 
present. The HBT data provide experimental evidence 
that favors a rapid breakup of QGP with a short time 
of hadron production 1211. and thus favors very fast, or 
sudden, hadronization [22, H^l- In this situation, a sim- 
ilar chemical non-equilibrium approach must be applied 
to the light quark abundance, introducing the light quark 
phase space occupancy 7^. This proposal made for the 
high energy SPS data helped improve the under- 

standing of RHIC200 hadron rapidity yield results 
and allowed a consistent interpretation of such data 
across the full energy range at SPS and RHIC200 [Tll |. 

Despite the abundant HBT evidence for fast hadroniza- 
tion, and the steady advances of lattice QCD _24-.27] pro- 
viding in past several years solid evidence of low QGP 
hadronization temperature, and further, despite our con- 
tinued effort showing that a good description of hadron 
abundances emerges with chemical non-equilibrium, the 
acceptance of light quark (m, d) chemical non-equilibrium, 
i.e., 7g 7^ 1, remained sparse. This situation can be re- 
lated to the misunderstanding of the meaning of lattice 
results by those insisting on light quark chemical equi- 
librium |13l - [l8l |. since these results predict hadronization 
and hadron production at a temperature above the lattice 
phase transformation temperature for 2-1-1 flavor. This 
is clearly impossible. 

The LHC2760 experimental environment has opened a 
new experimental opportunity to investigate in detail the 
SHM hadron production model. Precise particle tracking 
near to interaction vertex in the LHC- ALICE removes the 
need for off-line corrections of weak interaction decays. 



and at the same time vertex tracking is enhancing the 
efficiency of track identification, increasing considerably 
the precision of particle yield measurement [1, . All 
data used in the present work was obtained in this way 
by the LHC-ALICE experiment for Pb-Pb collisions at 
yj/VTV = 2.76 TeV, limited to the central unit of rapid- 
ity interval —0.5 < y < 0.5. The experimental particle 
yield results are reported in different collision central- 
ity bins according to the geometric overlap of colliding 
nuclei, with the 'smallest', e.g., 0-5% centrality bin cor- 
responding to the nearly fully ove rlap ping collision geom- 
etry. Geometric collision model [29] relates the central- 
ity trigger to the number of participating nucleons A^part 
which we use as our preferred centrality variable in what 
follows. 



Section |TT] presents our general method and approach 
to the particle multiplicity data analysis. Following a 
brief summary of the SHM methods in subsection III Al 
we describe in subsection IIIBI our centrality study of 
particle production based on the following data: for the 
0-20% centrality bin, we obtain the preliminary data 
from [2^. For the centrality study of particle pro- 
duction, we present in table Utt^'^, and p^ final yields 
as p resented in [s^l. The preliminary ratio (f>/K, is from 
|3l| ; these 7 data fields are binned in the same centrality 
bins and are used as presented. However, several other 
particle types require rebinning with interpolation and, 
at times, extrapolation which is further discussed in ap- 
pendix 1^ The (preliminary) data input into this rebin- 
ning for K*°/K~ and for 2A/(7r~-|- 7r+) are also taken 
from Ref. [31|. Using the preliminary enhancement fac- 
tors of S~, S"*", il~ , fl^ shown in Ref. [13, H^l, com- 
bined with yields of these particles for p-p-reactions at 
^/snn = 7 TeV as presented in Ref. [3J], we obtain the 
required yield input, see appendix 1X1 In subsection III C) 
we present both the particles fitted and predicted by 
SHM, including anti-matter clusters. 

In section IIIII we discuss the key physics outcome of 
the fits, i.e., the resulting SHM parameters as a function 
of centrality. We compare to the equilibrium approach in 
subsection nil Al We discuss the differences seen between 
the SHM variants and compare our results to our analy- 
sis of Au-Au collisions at y^SNN = 62.4 GeV at RHIC62, 
as it is a system we analyzed in detail recently [l3|. We 
obtain the bulk physical properties: energy density, en- 
tropy density, and pressure, as a function of centrality in 
subsection nil CI where we also address strangeness and 
entropy yields. This study is made possible since all SHM 
parameters are determined with minimal error in consid- 
eration of the precise experimental particle multiplicity 
result. We discuss how our results relate to the lattice- 
QCD study of QGP properties in subsection IIIIDI We 
close our paper with a short summary and discussion of 
all results in section ITVl 
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II. SHM AND PARTICLE PRODUCTION 

A. Generalities 

We use here SHM implementation within the SHARE 
program The SHM describes the yields of particles 
given the chemical freeze-out temperature T and over- 
all normalization dV/dy (as the experimental data are 
available as dN/dy). We account for the small asym- 
metry between particles and anti-particles by fugacity 
factors Xg, Xg and the light quark asymmetry A73, see 
Ref.Q- We further note that it is not uncommon to 
present the particle-anti-particle asymmetry employing 
the baryo-chemical and strangeness chemical potentials 
defined by 



^B=3TlnA, and fis ^ Tln{Xg/ Xg), 



(1) 



the 'inverse' definition of fig with reference to has 
historical origin and is source of frequent error. 

For each value of Xq, strangeness fugacity As is eval- 
uated by imposing the strangeness conservation require- 
ment (s) — (s) ~ 0. From now on, we omit the bra-kets 
indicating grand canonical average of the corresponding 
summed particle yield. The isospin fugacity factor A/3 
is constrained by imposing the charge per baryon ratio 
present in the initial nuclear matter state at initial in- 
stant of the collision. We achieve this objective by fitting 
these conservation laws along with particle yield data, us- 
ing the following form: 



s — s 



0.00 ±0.01, 



Q-Q 

- — = = 0.38 ±0.02. 

B-B 



(2) 

(3) 
(4) 



We believe that implementing conservation laws as data 
points with errors accounts for possibility that particles 
escape asymmetrically from the acceptance domain. 

In the LHC2760 energy regime, there is near symme- 
try of particle and anti-particle sector thus the chemi- 
cal potentials are hard to quantify. Therefore, the two 
constraints Eq.[2] and Eq.[3] alone were not sufficient to 
achieve smooth behavior of the chemical potentials as 
a function of centrality. We therefore impose as a fur- 
ther constraint a constant baryon number stopping per 
participating nucleon in the mid-rapidity region in the 
following form: 



b-b 



0.0054 ± 1%. 



(5) 



part 



We selected condition Eq.[S] since this was the variable 
which emerged in unconstrained fits as being most con- 
sistent. The value we selected is our estimate based on 
convergence without constraint to this value at several 
centralities. The alternative to this approach would have 



been to take across centrality a constant value of fiB- 
While this produces a good enough fit as well, this ap- 
proach was poorly motivated: the unconstrained fit re- 
sults produced rather random looking distribution of 
across centrality and thus did not present any evidence 
pointing towards a specific choice for ^g. While the ac- 
tual method of fixing matter-anti-matter asymmetry is 
extraneous to the main thrust of this paper, the value of 
fiB is of some relevance when considering predictions for 
anti-nuclei which we present further below. 

Our considerations include the already described phase 
space occupancy parameters 7s and 7^, where the light 
quarks q = u,d are not distinguished. We do not study 
7c here, in other words, we do not include in present 
discussion the charm degree of freedom. We note that 
there is no current experimental p_L-integrated charmed 
hadron yield information available. The integration of 
the phase space distribution is not yet possible due to 
uncertain low transverse momentum yields. 

Thus, in LHC2760 energy domain, we have at most 
4 = 7 — 3 independent statistical model parameters- 
constraints: seven parameters dV/dy,T,Xq,Xs,Xi3,jq 
and 7s constrained by the three conditions, Eq.[21 Eq.O 
and Eq.[5l to describe within SHM approach many very 
precise data points spanning in yield across centrality 
more than 5 orders of magnitude. We will show for com- 
parison results obtained setting arbitrarily 7^ = 1 (chem- 
ical semi-equilibrium fit, comprising 6 — 3 parameters- 
constraints) and than 7^ = 7s = 1 (chemical equilibrium 
fit, 5 — 3 parameters-constraints). 

Absolute yields of hadrons are proportional to one 
power of 7q for each constituent light quark (or anti- 
quark) and one power of 7s for each strange quark (or 
anti-quark). For example, 7^ enters non-strange baryon 
to meson ratios in the following manner: 



baryon(ggg) 7q 
meson{qq) 7^ 



meson y j 



(6) 



where q stands for either m or d quark and F is the inte- 
gral over all particle momenta of the phase space distribu- 
tion at freeze-out temperature: we always use exact form 
of relativistic phase space integrals. For strange hadrons, 
we must replace 7^ by 7s for each constituent s (and/or 
s) quark. Experimentally measured light baryon to me- 
son ratios (such as p/tt) strongly depend on the value of 
7g in a fit. Similarly, A{qqs)/7r{qq) cx 7s is very sensitive 
to the value of 7s . 

The value of 7^ is bound by appearance of a pion con- 
densate which corresponds to a singularity in the pion 
Bose-Einstein distribution function reached at the con- 
dition 



7o 



exp 



2T 



(7) 



This numerically works out for T = 138-160 MeV to be 
in range 7"'' = 1.63-1.525. On the other hand, there 
is a much more lax limit on the range of 7s , strangeness 
can increase very far before a particle condensation phe- 
nomenon limit is reached. 
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FIG. 2: (color online) Total ^-s a function of central- 
ity, as indicated in figure, for the total equilibrium (7, — 
7s = 1, ndf = 11), for the semi-equilibrium (7, = 1, 7^ / 
1, ndf = 10), and for chemical non-equilibrium (7, 7^ 1, 7^ 7^ 
1, ndf — 9) SHM. Open symbols represent the total for 
the 0-20% centrality bin depicted on left in figure [T] The num- 
ber of degrees of freedom for the three cases, respectively, is 
ndf = 11,10,9. 



B. Centrality study 

The input hadron yield data used in the fit to the 0- 
20% geometric centrality bin is shown in the 2nd column 
of table m The fit to this data set for the case of chemical 
equilibrium, where one forces 7^ = 7g — 1, was done 
in Ref. [l^ choosing a fixed value fiB = ^ MeV. In a 
first step, we compare to these results and follow this 
approach. However, we consider it necessary to apply 
strangeness and charge per baryon conservation by fitting 
Eq.[2]and Eq.[3]as two additional data points determining 
the corresponding values of chemical parameters /ig, fiis, 
a procedure omitted in the report Ref. ^2^ where fj-s = 
= was set. Naturally, the effect of this improvement 
is minimal, but it assures physical consistency. We show 
the values of Xtotai figured! see the large open symbols. 
The wider range of A'^part corresponding to the centrality 
bin 0-20% is shown in figure[5]as a horizontal uncertainty 
bars. 

In our detailed centrality dependent analysis, we use 
data in nine finer centrality bins, which we show in the 
third to eleventh and last column of table H] The bins 
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FIG. 3: (color online) Experimental charged particle yield 
pseudo-rapidity density dN/drj (blue squares), and geometric 
model relating to charged particle rapidity density dN/dy of 
only primary particles (violet circles) and including feed from 
strongly decaying resonances (red triangles) per participant 
pairs iVpart/2. 



are classified according to the average number of partic- 
ipants iVpart as a measure of centrality. This is a model 
value originating in the experimentally measur ed p seudo- 
rapidity density of charged particles dNch/drj [23], which 
we state in the third row of table HI We consider the 
consistency in figure [3] the experimentally measured 
dNch/dr] in the relevant participant bins [s^l is shown by 
square symbols, as well as our SHM results for rapidity 
density of charged particles dNch/dy emerging directly 
from QGP (i.e., primary charged hadrons) and the final 
yield of charged hadrons, as triangles, fed by the decay of 
hadronic resonances. In all cases, we show the yield per 
pair of interacting nucleons using the model value iVpart- 
While the primary charged hadron yield is well below 
the pseudo rapidity density of charged hadrons, the final 
yield is well above as is found in models accounting for 
the missing particle ID 1]. 

In figure [21 we see that about 50% of charged hadronic 
particles are produced by strong decays of heavier res- 
onances. We show, in figure [4l the ratio of primarily 
produced yield to the total yield for different particle 
species in the expected range of hadronization temper- 
atures. The dominant fraction, almost 80%, of tt and 
p yield originates from decaying resonances. This result 
demonstrates the difficulty that one encounters in the in- 
terpretation of transverse momentum spectra which must 
account for the decays and is thus, in a profound way, im- 
pacted by collective ffow properties of many much heav- 
ier hadrons [ll-Q. Conversely, this means that one can 
perform a convincing analysis of transverse momentum 
distribution only for hadrons, which do not have a signif- 
icant feed from resonance decays, such as or 0. This 
finding is the reason why we study the p± integrated 
yields of hadrons in exploration of the physics of the fire- 
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TABLE I: Table of data points we use as input for SHM fits; header of the table defines the centrality bins in three different 
ways. Hadron yields (dN/dy) and ratios at mid-rapidity \y\ < 0.5 for different centralities. Centrality as a function of A'part 
taken over from [29l ]. Errors are combined systematic and statistical errors added in quadratures where systematic errors are 
in general dominant and statistical errors are negligible. Values in brackets are interpolated data. See appendix [X] for details 
about data sources and how data are rebinned. 



Centrality 


0-20% 


0-5% 


5-10% 


10-20% 


20-30% 


30-40% 


40-50% 


50-60% 


60-70% 


70-80% 


{iVp„t> 


308 


382.8 


329.7 


260.5 


186.4 


128.9 


85.0 


52.8 


30.0 


15.8 


^ AT 1 

"-^ *ch / 




16U1 ± DO 


izy4 it 4y 


966 ± J ; 


649 ± z3 


426 lb 15 


261 ± 9 


149 ± 6 


76 ± 4 


35 ± 2 


7r+ 


562 ± 36 


733 ± 54 


606 ± 42 


455 ± 31 


307 ± 20 


201 ± 13 


124 ± 8 


71 ± 5 


37 ± 2 


17.1 ± 1.1 


7r~ 


560 ± 34 


732 ± 52 


604 ± 42 


453 ± 31 


306 ± 20 


200 ± 13 


123 ± 8 


71 ± 4 


37 ± 2 


17.0 ± 1.1 


K+ 


84 ± 5.4 


109 ± 9 


91 ± 7 


68 ± 5 


46 ± 4 


30 ± 2 


18.3 ± 1.4 


10.2 ± 0.8 


5.1 ± 0.4 


2.3 ± 0.2 


K- 


84 ±5.7 


109 ± 9 


90 ± 8 


68 ± 6 


46 ± 4 


30 ± 2 


18.1 ± 1.5 


10.2 ± 0.8 


5.1 ± 0.4 


2.3 ± 0.2 


K*" 


17.3 ± 4.2 




















K^'/KIO^ 




(188 ± 98) 


(196 ± 77) 


(209 ± 54) 


(227 ± 59) 


(247 ± 64) 


(269 ± 70) 


(295 ± 77) 


(326 ± 85) 


(361 ± 94) 


P 


25.9 ± 1.6 


34 ± 3 


28 ± 2 


21.0 ± 1.7 


14.4 ± 1.2 


9.6 ± 0.8 


6.1 ± 0.5 


3.6 ± 0.3 


1.9 ± 0.2 


0.90 ± 0.08 


P 


26.0 ± 1.8 
9.6 ± 1.4 


33 ± 3 


28 ± 2 


21.1 ± 1.8 


14.5 ± 1.2 


9.7 ± 0.8 


6.2 ± 0.5 


3.7 ±0.3 


2.0 ± 0.2 


0.93 ± 0.09 


4> 




109 ± 20 


116 ± 18 


117 ± 17 


128 ± 19 


120 ± 20 


123 ± 18 


123 ± 19 


119 ± 18 


119 ± 21 


A 


19.3 ± 2.0 




















A/tt 10^ 




(33.3±3.9) 


(34.2±4.0) 


35.3±4.1 


(36.4±4.3) 


(37.0±4.3) 


(37.1±4.4) 


(36.8±4.3) 


(36.0±4.2) 


(34.7±4.1) 


E- lO--' 


3230 ± 350 


(3970 ± 440) 


(3370 ± 370) 


(2580 ± 280) 


(1760 ± 190) 


(1160 ± 130) 


(716 ± 79) 


(407 ± 45) 


(196 ± 22) 


(75 ± 8) 




3040 ± 330 


(3820 ± 420) 


(3270 ± 360) 


(2530 ± 280) 


(1760 ± 190) 


(1180 ± 130) 


(737 ± 81) 


(423 ± 47) 


(203 ± 22) 


(71 ± 8) 


SI" 10^ 


570 ± 100 


(780 ± 188) 


(620 ± 155) 


(450 ± 110) 


(290 ± 70) 


(180 ±40) 


(102 ± 26) 


(54 ± 14) 


(25 ± 6) 


(10 ± 3) 


n+ 10^ 


580 ± 110 


(760 ± 184) 


(610 ± 152) 


(450 ± 110) 


(290 ± 70) 


(180 ±40) 


(102 ± 26) 


(55 ± 14) 


(25 ± 6) 


(10 ± 3) 
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FIG. 4: (color online) Fraction of primary hadrons produced, 
normalized by their final yield, which consists of primary pro- 
duced hadrons and the feed from strong decaying resonances, 
for particles indicated on the right margin, as function of 
hadronization temperature in single-freeze-out model. 



ball particle source. Moreover, we believe that 'blast- 
wave' model fits to p± hadron spectra are only meaning- 
ful for the 51 or hadrons. 

The centrality binning, which differs for different par- 
ticles considered, requires us to use several interpolated 
and even some slightly extrapolated experimental results, 
which procedure we discuss in depth below, and in ap- 
pendix [Xj In our fits, we choose to use the centrality bins 
with the largest number of directly determined experi- 
mental data, minimizing the potential error originating 
in our multi-point interpolation. A few particles appear 
more than once in our data set (as a yield and/or in a 
ratio). However, to prevent duplicity, we always fit every 
particle measured just once. 



In order to show that the finer centrality binning mat- 
ters, we have already shown the 5-10% centrality bin 
(which contains only close extrapolation and interpola- 
tion, see open symbols) in the right panel of figure [T] 
The fit has the same set of particles as the 0-20% cen- 
trality bin seen on left in figure [TJ however, some of these 
particles enter the finer binned fit in ratios. If the out- 
come of the fit as a function of centrality is even a small 
variation in fitted parameters (other than normalization, 
i.e., volume), we expect and we find that the 5-10% cen- 
trality bin, which describes a much smaller participant 
.^part range, leads to smaller compared to the wide 
0-20% case. However, the stability of the fit parameters 
implies that much of the improvement is attributable to 
the revision in the input data set, see appendix El Again, 
the non-equilibrium SHM variant shows its superiority by 
providing a factor of 4 smaller x^- 

We perform a fit to the entire data set with all three 
SHM approaches and compare the resulting as a func- 
tion of iVpart in figure [2J The solid squares represent the 
chemical non-equilibrium SHM (7^ ^ 1,7^ 7^ 1), the solid 
circles represent the semi- equilibrium SHM (7, = 1 , 7^ 7^ 
1) and solid triangles represent the full equilibrium SHM 
(79 = 7s = !)• The range of centrality is indicated by the 
horizontal bars. Considering most central bins, we note 
in figure [5] that allowing 7^ 7^ 1 can reduce the total 
of the fit by more than a factor of 3 compared to semi- 
equilibrium, and more than a factor of 5 comparing full 
non-equilibrium with full equilibrium. 

As a last step, we verify if there is a special value of 
the parameter 7^ of particular importance. To this end, 
we have evaluated the x^/ndf of the fit as function of a 
given fixed 7^ within a range 7^ G (0.95,7"'*). This 
profile curves, seen in figure O all pass 7q — 1 smoothly, 
therefore 7g = 1 has no special importance for the SHM. 
However, fits to data in all centralities decrease in x^ as 
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FIG. 5: (color online) x^/ndf profile as a function of 7q for 
all studied centralities. 



7^ increases, they all point to best fit value of 7g near 
the critical value of Bose-Einstein condensation given by 
Eq.[71 We have to note that the most peripheral bin 
(70-80%, iVpart = 15.7) analyzed here shows consider- 
ably lower for equilibrium and semi-equilibrium. This 
could indicate that a small fireball created in peripheral 
heavy-ion collisions by little more than a dozen partici- 
pants may not always enter the QGP phase. 
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C. Particle yields 

We compare input to the resulting particle yields 
graphically in figure [51 We fit 13 particles, counting 
anti-particles, which in the figure cannot be visually dis- 
tinguished as an independent input or output data, and 
ratios A/tt, K*/K and (j)/K. For these ratios, the relevant 
yield outputs A, K* and </> are shown. Direct comparison 
of the input A/tt, K*/K and t/f/K ratios to the output is 
presented in figure [71 note that (j)/K ratio is available as 
experimental data point in all centrality bins. The fitted 
output yields are stated also in the top portion of tablellll 
and ratios are given just below allowing for comparison 
with the input values. 

Our fit results appear as open circles in figure [HI at 
times completely overlaying the input data, full symbols. 
For the A, the dotted line guides the eye, since the ac- 
tual fit is to the ratio A/tt shown in figure [71 no absolute 
A data are available, in absence of absolute yields, only 
open circle, i.e., the fitted value, is shown in figure[6l Sim- 
ilar situation arises with (j) and K*, where data are not 
available, but we fit (j)/K and K''*/K. One can see that 
SHM generated results follow closely both the experimen- 
tal data available, and the interpolation dashed lines for 
each particle, and that each interpolation curve passes 
through the experimental data points shown in full sym- 
bols, or at worse, the error bars if these are larger than 
the symbol. 

Even so, we note in figure [51 that our interpolation 



FIG. 6: (color online) Full symbols are the experimental data 
points. Open symbols represent the outcome for our chemical 
non-equilibrium SHM fit to LHC2760 as a function of central- 
ity, i.e., A^part- Dashed lines are outcome of our data inter- 
polation of experimental yields available (see appendix . 
Dotted lines connecting A, K* and (p SHM output values are 
presented to guide the eye; since ratios of these particles were 
used in our fit, the data are shown in figure [7| 



for H. shows a slightly different systematic shape (dashed 
line) compared the fit results (open symbols) or the be- 
havior of the other particles. In other words, we see that 
other particles 'predict' a yield of that 'looks' more 
like that of other particles, while the four data points 
lead to a centrality distribution that is slightly different. 
More precise data will without doubt offer a resolution 
to this slight tension in our interpolation/extrapolation. 
The hadron yields we find are also stated in table [ill 
Aside from the yields, we show there frequently quoted 
ratios of particle yields, e.g., we find p/7r+ ~ 0.05. We 
will return to discuss this ratio in subsection IIII Bl 

Figure[7lhas the largest differences between theory and 
experiment. In case of the A/tt ratio, we see a systematic 
within error bar under-prediction at all centralities. For 
K*/K, we see within the error bar a different slope of the 
fit as a function of centrality. The question can be asked 
if these differences of fit and results indicate some not 
yet understood physics contents. However, we are within 
error bars and such data-fit difference must be expected 
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TABLE II: Table of hadron yield output; header of the table defines the centrality bins in three different ways. Top section 
of the table shows fitted yields dN/dy of hadrons entering the fit at LHC2760 obtained in the chemical non-equilibrium SHM. 
Next below are three ratios that are actually included in the fit (rather than the yields of A, K"*, </!)), followed by the ratios of 
hadron yields that can be formed from the stated results, stated here for convenience of the reader. In the two lower sections 
of the table, there are predicted yields of yet unmeasured hadrons and at the very bottom, we show predicted yields of light 
anti-nuclei scaled up by factor 1000 (and by 10® for anti-Helium). Note that yield of matter particles is nearly the same. 
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30-40% 


40-50% 
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382.8 


329.7 


260.5 


186.4 
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30.0 
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696 
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35.8 
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0.953 




20.8 


27.0 


22.6 


17.0 


11.7 


7.56 


4.68 


2.67 


1.36 


0.587 


4> 


11.2 


14.4 


12.1 


9.12 


6.23 


3.99 


2.40 


1.33 


0.632 


0.255 


A 


17 9 




lo.u 


14 9 


Q 7'i 


u.^o 


A nn 


9 '^1 


1 IS 
1 . lo 


n ^9n 




o.uo 


o.ou 


% 9"^ 


9 AQ 

Z.4;cy 


1 7n 

1 . / u 


111 
1.11 




n '^77 


n 1 81 

U. lol 


U.U t 10 




3.00 


3.80 


3.22 


2.43 


1.68 


1.08 


0.664 


0.371 


0.179 


O OTOi? 

U.Uyzd 


n~ 10^ 


529 


663 


561 


435 


297 


193 


115 


62.6 


28.5 


11.0 
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32.9 
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fi/TT 10^ 


1.007 


0.952 


0.976 


1.010 


1.019 


1.019 


0.973 


0.920 


0.795 


0.671 


K°77r 10^ 


39.6 


38.8 


39.3 


39.5 


40.0 


39.9 


39.7 


39.2 


37.8 


35.9 


p/K 


0.300 


0.303 


0.302 


0.306 


0.309 


0.322 


0.336 


0.354 


0.385 


0.426 


■n 


61.0 


79.7 


66.3 


49.8 


33.9 


21.9 


13.4 


7.64 


3.90 


1.72 


p(770)° 


38.9 


51.3 


42.5 


32.1 


21.9 


14.5 


9.14 


5.35 


2.87 


1.34 


a;(782)0 


35.1 


46.4 


38.4 


29.0 


19.8 


13.1 


8.28 


4.85 


2.61 


1.22 


A(1232)++ 


4.98 


6.57 


5.46 


4.15 


2.86 


1.92 


1.23 


0.734 


0.402 


0.191 


E*(1385)- 


2.08 


2.70 


2.26 


1.72 


1.18 


0.785 


0.492 


0.284 


0.146 


0.065 


A*(1520) 


1.09 


1.41 


1.18 


0.907 


0.625 


0.418 


0.264 


0.153 


0.0795 


0.0355 


H*(1530)- 


1.02 


1.30 


1.09 


0.84 


0.58 


0.378 


0.230 


0.129 


0.0626 


0.0258 


^HIO^ 


74.7 


98.1 


81.9 


62.6 


43.6 


29.4 


19.5 


11.8 


6.53 


3.16 


iHio^ 


0.478 


0.601 


0.506 


0.397 


0.279 


0.191 


0.128 


0.0773 


0.0415 


0.0193 


^HIO^* 


1.64 


2.13 


1.79 


1.39 


0.983 


0.677 


0.468 


0.290 


0.166 


0.083 


^HilO^ 


1.64 


2.14 


1.79 


1.39 


0.986 


0.679 


0.469 


0.291 


0.166 


0.083 


■'HelO® 


5.87 


7.57 


6.41 


5.04 


3.64 


2.56 


1.85 


1.18 


0.697 


0.362 



and is allowed given a large data sample and potential for 
experimental refinement of these two preliminary data 
sets involving K* and A. We recall that at RHIC200, 
K*/K ratio was flatter, more independent of centrality 
with a value close to the most central LHC2760 result. 
We also note that we did not yet study how the charmed 
hadron decay particles influence the fit. 

Predictions for the six hadron yields 77, p", A(1232)++, 
A*(1520), S]*(1385)-, S*(1530)0, are shown in figure [S] 
as a function of centrality, these results are stated in the 
lower portion of table HT] We further show five differ- 
ent species of (strange) anti-matter, from anti-deuteron 
to anti-alpha, including anti-hypertriton, appropriately 



scaled to fit into display of figure HI Our predictions 
of these composite objects should serve as a lower limit 
of their production rates: fiuctuation in the QGP ho- 
mogeneity at hadronization, and recombinant formation 
after hadronization may add contributions to the small 
SHM yield, see here the corresponding RHIC result [35| . 



9 



0.40 
0.30 

0.20 
0.15 

0.10 
0.07 
0.05 

0.03 



1 K7K 

If* f 1 I 
^ (^/K 



T A/ti 



50 100 150 200 250 300 350 400 
Npart 

FIG. 7: (color online) Full symbols: experimental data with 
errors for K*/K, (^/K, and A/tt as a function of centrality, i.e., 
number of participants A'part- Lines are the interpolation, re- 
spectively, extrapolation except for the case (/I>/K, where dot- 
ted line guides the eye. Open circles represent the resulting 
fit value for each ratio. The (blue) shaded band shows the 
error input used obtaining the interpolated values of K*/K 
ratio. 



III. PARTICLE SOURCE AND ITS 
PROPERTIES 

A. Statistical parameters 

In figure [51 we depict the LIIC2760 statistical param- 
eters as a function of collision centrality and compare 
these LHC2760 results with those we have obtained at 
RHIC62 shown by open symbols. In all three pan- 
els of figurejni we show parameter errors evaluated by 
SHAREv2 J4| employing the MINOS minimization rou- 
tine. One can see that the parameter values for chemical 
non-equilibrium are defined better than for the case with 

We present LHC2760 hadronization parameters for the 
non-equilibrium SHM case also in the top section of ta- 
ble mil In the top frame figure [5^, we see the particle 
source volume dV/dy, in the middle frame figure ^}p, 
the chemical freeze-out temperature T, and in the bot- 
tom frame figure [SJ;, the phase space occupancies — 
the different variants are distinguished by superscripts 
'neq' (non-equilibrium, that is 7^ 7^ 1,7s 7^ 1) , 'seq' 
(semi-equilibrium, 7^ = 1,7^ 1). and 'eq' (equilib- 
rium, 7q = 1,7s = 1). To compare with the semi- 
equilibrium SHM variant, wc show the ratio 7s"°'^'' /t?""'*'' i 
a ratio which helps to quantify the strangeness to light 
quark enhancement. This is to be directly compared with 
the semi-equilibrium strangeness phase space occupancy 
7s^°^\ given fixed 7^'''"''' = 1. 

For the LHC2760 data, the SHM forcing chemical equi- 
librium of light quarks (i.e., 7^ = 1 with either 7^ = 1 or 
7s 7^ 1) have a very similar volume dV/dy, and similar 
chemical freeze-out T as shown in figure [S^jb, respec- 
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FIG. 8: (color online) Particles predicted by the chemical 
non-equilibrium SHM fits. Predictions for hadron yields (full 
symbols) are above anti-nuclei yield predictions (open sym- 
bols), which have been multiplied by suitable factors. Lines 
guide the eye, points are actual predictions for each of the 
nine centralities we analyzed. 



tively, with nearly overlapping lines for dV/dy. In the 
non-equilibrium approach, dV/dy is reduced by about 
20-25%, and the freeze-out temperature T by 10% com- 
pared to the equilibrium SHM variant. Compared to the 
RHIC62 results ^ (open symbols) the LHC2760 vol- 
ume dV/dy is up to a factor 4 larger while the LHC 
hadronization temperature T is 2-5 MeV lower. Thus 
the main origin for the much larger particle multiplicity 
found within the SHM interpretation of hadron produc- 
tion data is the considerably increased transverse dimen- 
sion of the fireball at time of hadronization. We under- 
stand this growth of volume as being due to a more vi- 
olent transverse fireball expansion, driven by the greater 
initial energy density formed in LHC2760 heavy ion colli- 
sion and corresponding initial pressure. This dynamical 
expansion leads to a greater supercooling, which explains 
the systematic decrease of the freeze-out temperature at 
LHC2760 compared to RHIC62. 

The freeze-out temperature T at LHC2760 decreases 
when considering more central collisions, see figure ^jp. 
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TABLE III: Top section shows chemical non-equihbrium SHM fit parameters dV/dy, T, 7,, 7^ and Xtotai with ndf (number 
data less number of parameters) obtained in each centrality bin. Bottom section presents fireball bulk properties in each bin; 
energy density e, pressure P, entropy density cr, strangeness per entropy content s/S, and entropy at LHC2760 compared to 
RHIC62, Slhc/Srhic. 



Centrality 


0-20% 


0-5% 


5-10% 


10-20% 


20-0% 


30-40% 


40-50% 


50-60% 


60-70% 


70-80% 


(Afpart) 


308 


382.8 


329.7 


260.5 


186.4 


128.9 


85.0 


52.8 


30.0 


15.8 


dV/dy [fm3] 


2463 ±6 


3304±469 


2715±81 


2003±47 


1337±173 


853.9±5.9 


512.2±70.1 


289.4±5.5 


149.8±5.0 


66.9±9.7 


T [MeV] 


138.3±0.0 


138.0±0.0 


138.1±0.0 


138.6±0.0 


139.0±1.7 


139.9±0.0 


140.7±1.8 


141.5±0.0 


142.3±0.5 


143.2±2.0 


7? 


1.63±0.00 


1.63±0.00 


1.63±0.08 


1.63±0.00 


1.62±0.08 


1.62±0.05 


1.62±0.08 


1.61±0.00 


i.eiito.oo 


1.60±0.09 


7s 


2.08±0.00 


2.05±0.17 


2.06±0.13 


2.06±0.13 


2.05±0.16 


2.00±0.01 


1.95±0.16 


1.88±0.05 


1.75±0.19 


1.62±0.14 


X?otal/ndf 


9.51/9 


3.87/9 


3.52/9 


3.94/9 


3.35/9 


4.73/9 


4.55/9 


5.65/9 


6.13/9 


4.09/9 


e [GeV/fm^] 


0.462 


0.453 


0.457 


0.467 


0.476 


0.487 


0.505 


0.516 


0.521 


0.527 


P[MeV/fm3] 


78.5 


77.1 


77.7 


79.1 


80.5 


82.3 


85.1 


86.8 


87.9 


89.2 


a [fm-3] 


3.20 


3.14 


3.17 


3.23 


3.28 


3.36 


3.46 


3.53 


3.56 


3.60 


s/S 


0.0299 


0.0295 


0.0297 


0.0297 


0.0298 


0.0294 


0.0290 


0.0284 


0.0272 


0.0257 


•S'lhc/Srhic 


3.07 


3.23 


3.10 


2.93 


2.75 


2.56 


2.33 


2.06 


1.74 


1.27 



This can be interpreted as being due to a deeper su- 
percooling of the most central and most energetic colli- 
sion systems. We can extrapolate the freeze-out tem- 
perature to A^part = in the figure to set an upper 
limit on hadronization temperature at LHC2760, Thad 
145 ± 4 MeV, applicable to a small (transverse size) fire- 
ball. This, then, is the expected hadronization temper- 
ature without supercooling. Excluding, in figure |91d, the 
most peripheral T-fit point for RHIC62, which does not 
have a good confidence level, we see that T at RHIC62 
converges towards the same maximum value as we found 
at LHC2760, thus confirming the determination of Thad 
as the common hadronization temperature without su- 
percooling. 

We show the phase space occupancies 7q,7s in fig- 
ure [9]:. We note that the LHC2760 fit produces nearly a 
constant 7^ as a function of centrality. However, 7^ (and 
respectively Js/^q) decrease for more peripheral collisions 
towards unity suggesting that these two flavors approach 
chemical equilibrium for system of small transverse size. 
A similar situation for peripheral collisions was observed 
for RHIC62. However, at RHIC62, we see a strong cen- 
trality dependence of 7s and hence 7s /7g- This rapid 
rise of the RHIC62 7s as function of centrality can be at- 
tributed to the buildup of strangeness in QGP formed at 
RHIC62, which is imaged in the later produced strange 
hadron yield. Note that, omitting the most peripheral 
RHIC62 point, the peripheral 7, is nearly the same as at 
LHC2760. The small difference can be attributed to the 
smaller allowed value of 7^ for the slightly higher value 
of T seen at RHIC62. 

We have executed all our fits allowing for the presence 
of the chemical potentials (Eq.[T]) characterizing the slight 
matter-anti-matter asymmetry present at LHC2760. The 
quality of the fit is not sufficiently improved including 
effectively one extra parameter (/is, since fis is fixed 
by strangeness conservation) to assure that the uncon- 
strained results for fiB sue convincing. As mentioned in 
section Hi A[ we smooth the centrality dependence of /is 



by introducing baryon stopping fraction at mid-rapidity, 
that is imposing Eq.[5] as an additional data point, a 
value that we saw a few times in the data without in- 
troducing this constraint. This constraint leads to the 
chemical potentials fiB and /is presented in figure [TUl 
with the baryochemical potential 1 < /^s < 2.3 MeV and 
fis — 0.0 ± 0.5 MeV for all centralities, values an order 
of magnitude smaller than at RHIC62 and RHIC200. As 
we can see, even with the constraint, there are two cen- 
tralities which do not agree with the trend set by the 
other seven data points. Data shown in figure [TUl are not 
defined well enough to argue that we see a decrease of 
baryochemical potential with increasing centrality, since 
this outcome could be result of the bias we introduced. 
However, we think that for the most central collisions at 
LHC2760 there is some indication that fiB — 1.5 MeV. 

B. p/tt ratio and chemical (non-)equilibrium 

The key difference between the three SHM approaches 
are the values of jq^s, as seen in figure[9};. In section Hi A[ 
we argued that the baryon to meson ratio, e.g., p/tt, is 
directly proportional to 7^ and this can be used to dis- 
tinguish between the three SHM approaches. This ratio 
is a big problem for the equilibrium SHM [§]. We wish 
now to quantify this result within our approach and to 
show that, within the chemical non-equilibrium SHM, the 
problem is solved. 

For this purpose, we redo all fits but making this ratio 
more explicit in the data analysis. Specifically, first we 
evaluate p/tt ratio based on the yields of p and tt seen in 
table U 

^ - (8) 

TT TT + 7r+ 

We estimate the error of the p/tt ratio by adopting the 
relative error of p/tt from that is 6.5%. We include 
this new data point, p/tt ratio, in the fit. Note that this 
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FIG. 9: (color online) SHM parameters as a function of cen- 
trality, i.e., number of participants A'^part, presented for the 
three different levels of chemical equilibrium, and compared 
to the chemical non-equilibrium SHM RHIC62 results (open 
symbols, dashed lines). All lines guide the eye. From top to 
bottom; (a) dV/dy — note that the volume for both equi- 
librium and semi-equilibrium SHM variants is so close that 
symbols overlap; (b) T, the chemical freeze-out temperature 
(semi-equilibrium symbols are offset in order to separate them 
from equilibrium); (c) phase space occupancies 7s"°''\ 75"°''' 
and for comparison with equilibrium 7!''°'''; we also present 
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increases the relative importance of p and tt compared to 
the other particles included in the fit. Open symbols in 
the top panel of figure [11] depict the data and full symbols 
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FIG. 10: (color online) Scatter plot of fitted chemical poten- 
tials. 



show the resulting output values obtained when we refit 
with enlarged data set that includes the p/tt ratio. There 
is a minimal change in statistical parameters and physical 
properties of the fireball which we do not restate. In 
bottom panel of figure [TTl we show Xtotai- 

Even with the increased importance of p/vr, the chemi- 
cal non-equilibrium SHM works very well. However, SHM 
with fixed 7^ = 1 have increased difficulties describing 
this ratio, that is there is systematic 1.5-2 s.d. differ- 
ence of the fit result and data and the value of Xtot^i 
large. When compared to the Xtotai obtained without the 
added p/tt in figure [21 the non-equilibrium variant shows 
nearly the same values of Xtotai ^'^^ centralities, the 
p/tt ratio is a natural outcome of the non-equilibrium 
approach. On the other hand, SHM approaches with 
7q = 1 show additional systematic increase in by a 
factor of ~ 1.3-1.5 for all centralities. This means that 
p/tt data are in conflict with the hypothesis 7^ 1. This 
demonstrates that the hypothesis of chemical equilibrium 
of light quarks is incompatible with the baryon to meson 
ratio at LHC2760 and 7^ ~ 1.6 is needed, in order to de- 
scribe the LHC data. This finding is in agreement with 
the RHIC200 data [s^, where the importance of the p/tt 
ratio was noted. 

To compare p/tt ratio with our predictions, recall that 
the picture of universal hadronization condition with uni- 
versal hadronization pressure P = 82 ± 5 MeV has been 
advanced by our group [ll|, [s^, HI]. For this favored 
hadronization condition, the p/tt ratio is predicted in ta- 
ble II of Ref. [39] to be p/tt = 0.047± 0.002, which agrees 
practically exactly with the experimental result shown 
in figure [TT] top frame. The ALICE collaboration [s^] 
presents the mechanism of post-hadronization proton- 
antiproton annihilation in order to justify the small p/tt 
ratio, as compared to result of equilibrium SHM models. 
Aside of annihilation, there are pp formation events. The 
significantly larger abundance (and therefore also den- 
sity) of heavy mesons compared to nucleons, see table [TT] 
implies that mesons can be an effective source of nucleon 
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FIG. 11: (color online) Top panel: data and SHM description 
of p/tt ratio fitted together with all other data within the 
three SHM approaches, as a function of centrality. Bottom 
panel: resulting Xtotai for Ihe three variants. See section UlIBI 
for details. 



pairs in reactions such as p + p i — > p + lu and many 
other relevant reactions, see table II in Ref. 40j. ALICE 
collaboration notes, that p/tt ratio modification after an- 
nihilation should disappear in most peripheral collisions 
due to smaller volume. We will now quantify this ef- 
fect showing how this fade-out of the annihilation effect 
would work as a function of centrality. We show that 
given the constant p/tt ratio in a wide range of centrali- 
ties, figure [m the effect of post-hadronization change of 
p/tt ratio must be negligible. 

In order to establish the centrality dependence of post- 
hadronization nucleon yield changing reactions, we eval- 
uate the total number of pp annihilation events. This 
number is obtained by integrating annihilation rate over 
history of the post-hadronization matter expansion 



A^annih = / dtNp{t) Pp{t)ae,nnihV , 



(9) 



where v is the relative velocity of p and p. The three- 
dimensional dilution of the density can be modeled as 



Ppit) 



Pv 



(1 



(10) 



{L)^[{dV/dy)/iAn/3)]'/\ 



(L) is the magnitude of the fireball size, and (uflow) — 
0.6-1 c is the velocity of the fireball expansion, in both 
cases averaged over the fireball complex three dimen- 
sional geometry. 



The initial density at time of hadronization is obtained 
from our hadronization study: 



Ppith) 



h h _ dNp/dy 
Pp- Pp 



dV/dy 



(11) 



In a wide range of low relative e nerg ies, which are relevant 
here, the event cross section is 



CTcvcnt = CTannihw/c ~ 46 mb. 



(12) 



Neglecting the depletion of nucleons (i.e., Np{t) ~ iV^), 
we find, combining Eq.[3] with Eg. llOl the ratio of anni- 
hilated (anti)protons to their total yield A^annih/^^ and 
proton mean path before it annihilates Lcvont: 



C 



(13) 



The upper three lines, in figure [121 show Lcvcnt for 
the three models of hadronization (equilibrium, semi- 
equilibrium and non-equilibrium) as a function of cen- 
trality. The colored band in figure [12] represents the er- 
ror originating from the freeze-out temperature T un- 
certainty (see figure [9]) . Note that the non-equilibrium 
model has much smaller parameter errors, so Lcvcnt is 
defined more precisely. Since the event reaction cross 
section for annihilation is well measured and nearly con- 
stant, it does not introduce any additional uncertainty 
to Lovont- The bottom three lines, in figure [Hj (semi- 
equilibrium and equilibrium lines overlap, since dV/dy is 
very similar in these two cases), show how the size {L) 
of the system changes as a function of centrality. Espe- 
cially for peripheral collisions, we see that (L) <C Lovont- 
The ratio of both length scales provides a measure of the 
fraction of protons that can be annihilated. 
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FIG. 12: (color online) Full symbols: antiproton annihilation 
event path Lcvcnt for the three SHM models as a function of 
centrality; open symbols: fireball size scale (L). Shaded error 
band represents error margin originating from the uncertainty 
of T (semi) equilibrium SHM. For readability, we omit much 
smaller error band for non- equilibrium SHM. 
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As seen in figure [T^l from central to semi-peripheral 
(-^part — 85) collisions, the ratio of both lengths nearly 
doubles. This means that the annihilation fraction drops 
in semi-peripheral collisions to about half of the most 
central value. However, the measured ratio p/tt is nearly 
constant over this range, increasing from 0.046± 0.003 to 
0.050 ±0.003. We interpret this as experimental evidence 
that the net effect of pp formation and annihilation is in- 
significant. Therefore, the annihilation of pp pairs cannot 
serve as the explanation of the disagreement between the 
equilibrium SHM and observed small value of p/tt ratio. 

Our estimate of the annihilation effect based on Ea.[T51 
and the result seen in figure [T^] is consistent with the 
annihilation effect reported in Ref. [1^, where detailed 
balance reactions forming pp were not considered. In 
this work, p/tt rises to p/tt = 0.058 already in the 20- 
30% centrality bin (iVpart = 185), which is more than 
3 s.d. above experimental data (see figure [11]). This 
shows that pure annihilation afterburners cannot explain 
a reduction of p/tt ratio from equilibrium SHM value. 
Therefore, we disagree with the remarks of the ALICE 
collaboration [soj , which, assuming the validity of chem- 
ical equilibrium, considers the pp annihilation process as 
a viable path to understand their experimental result for 
the p/tt ratio. 



C. Fireball bulk properties 

In order to obtain the bulk physical properties of the 
source of hadronic particles, we use exactly the same set 
of particles and the same assumptions about their prop- 
erties as we employed in the fit procedure. Therefore, the 
physical properties we determine are consistent with the 
particle yields that originated our fit. In other words, we 
sum the energy, entropy, etc. carried away by the ob- 
served particles, adding to this also what the predicted, 
but unobserved particles, 'contain'. 

The bulk physical properties of the hadronizing fire- 
ball, that is energy, pressure, entropy and strangeness 
per entropy content are shown in the bottom part of ta- 
ble mil and in figure [T3] where shaded domains show our 
error estimate. Solid symbols are results of the fit, lines 
guide the eye. In our SHAREv2 fit with MINOS mini- 
mization, the largest uncertainty seen in table IIIII is the 
7s and dV/dy error, see figure [H other statistical bulk 
properties have relatively insignificant errors. All fit er- 
rors propagate into the properties seen in figure [T^ Since 
in figure [H] we consider densities, the error in volume does 
not affect these values. Therefore, by recomputing the 
properties of the fireball shifting alone the value of 7s 
within one s.d., we obtain a good error evaluation in the 
measurement of the bulk physical properties shown in fig- 
ure [131 The point that stands out with very small error 
is at iVpart — 130. This anomaly is due to accidental 
appearance of a sharp minimum in the highly non-trivial 
7-dimensional parameter space. 

We are interested to study the bulk properties of the 
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FIG. 13: (color online) Bulk properties of the fireball as a 
function of centrality: from top to bottom, energy density e 
(purple triangles), the entropy density a (red) circles, scaled 
down by factor 10, and the hadronization pressure 3P (blue) 
squares, LHC2760 values shown with full symbols, RHIC62 
with open symbols for comparison. Shaded areas show our 
estimate of systematic error arising from the uncertainty of 
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source of hadrons in order to test the hypothesis that a 
QGP fireball was the source of particles observed. For 
this to be true, we must find appropriate magnitude of 
bulk properties consistent with lattice results, and at the 
same time, a variation as a function of centrality that 
makes good sense. We observe in figure [T3| a smooth and 
slow decrease of energy density e (top) , entropy density a 
(middle) and hadronization particle pressure P (bottom) 
as a function of centrality. This slow systematic decrease 
of all three quantities is noted in particular comparing 
to RHIC62 (open symbols) , where the properties seem to 
vary less. This maybe interpreted as an effect of volume 
expansion at LHC leading to larger supercooling for larger 
systems. 

The local thermal energy density of the bulk is the 
source of all particles excluding the expansion flow kinetic 
energy. The value we find is £ ~ 0.50 ± 0.05GeV/fm'^ 
in the entire centrality range. Nearly the same value is 
found within the chemical non-equilibrium approach for 
RHIC62 [il] and RHIC200 We note that e assumes 
the smallest value for the most central collisions, see ta- 
ble IIIII and figure [T31 The hadronization pressure P and 
entropy density a are also decreasing for more central 
collisions, which is consistent with our reaction picture 
of expanding and supercooling fireball — the larger sys- 
tem supercools a bit more. The error band is (as for e) 
based on 7s uncertainty. 
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FIG. 14: (color online) Entropy yield dS/dy at LHC2760 and 
at RHIC62 as a function of centrality showing power law fit 
parameters in the legend. Colored band represent uncertainty 
based on 7s fit uncertainty. Error bars arise from error in the 
volume dV/dy. 



In the last row of table IIIIl we show that entropy yield 
at LHC2760 is more than 3 times greater than obtained 
at RHIC62. The entropy yield dS/dy, as a function of 
participant number, is shown in figure lT^ and the notable 
feature is that the power law parametrization displays a 
nearly linear dependence at RHIC62 while at LHC2760 a 
strong additional entropy yield, associated with the faster 
than linear increase, is seen: dS/dy oc N^;^^^'^. Most of 
the entropy is produced in an initial state mechanism 
which remains to be understood and our finding of the 
nonlinear entropy growth with iVpart adds to the entropy 
production riddle and important observational result. 

However, at LHC2760, one expects a component in 
the entropy count arising from the inclusion of the de- 
cay products of heavy charmed hadrons in the hadron 
yield. This entropy component is different from entropy 
produced in initial reactions, this is the entropy aris- 
ing from hard parton collision production of charm, and 
post-hadronization decay of charmed hadrons. It is un- 
likely that the non-linearity of the entropy yield is due to 
this phenomenon as one can easily see that the required 
charm yield would be very large. We will return, in near 
future, to this question. The uncertainty of entropy de- 
picted in figure [14] as a shaded band is based alone on 7s 
variation, as was obtained for other physical properties 
in figure 1131 A further error due to variance in dV/ dy is 
shown as a separate error bar. Where it is invisible in 
the LHC2760, it is hidden in symbol size. 

We turn now to study strangeness per entropy s/S = 
{ds / dy) / {dS / dy) in the source firebaU. We are inter- 
ested in this quantity since both entropy and strangeness 
yields are preserved in the hadronization process. There- 
fore, by measuring s/S, we measure the ratio of strange 
quark abundance to total quark and gluon abundance 
which determines the source entropy, with a well known 
proportionality factor. For the presently accepted small 
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FIG. 15: (color online) Top frame: strangeness per entropy 
s/S content of the fireball at LHC2760 (filled squares) and 
at RHIC62 (open squares) as a function of centrality; bot- 
tom frame: the thermal energy cost to make a strange- 
anti-strange quark pair. Colored bands represent uncertainty 
based on 73 uncertainty. 



strange quark mass ms(/i — 2GeV) = 95 ± 5 MeV |43| . 
the predicted value shown in figure 5 of Ref. [4J| is 
s/S ~ 0.0305 ± 0.0005. Finding this result in our LHC 
data analysis is necessary in order to maintain the claim 
that the source of hadrons is a rapidly disintegrating 
chemically equilibrated QGP fireball. 

In the top panel of figure 1151 we show the strangeness 
per entropy s/S in the source fireball. The solid squares 
are for the LHC2760, and open symbols for RHIC62. We 
see that s/S saturates at s/S ~ 0.030 at LHC2760, 
a value reached already for A^part > 150, thus for a 
smaller number of participating nucleons than we found 
at RHIC62, and which value remains constant up to the 
maximum available iVpart- This agrees with equilibrated 
QGP hypothesis and suggests that the source of hadrons 
was in same conditions for a wide range of centrality. 

This constant s/S value as function of centrality can 
be interpreted as evidence of chemical equilibrium for 
a QGP source: the strangeness yield normalized to all 
quark and gluon yield inherent in S can be constant only 
if dynamical processes find a chemical balance for the 
differently sized fireballs. The value s/S = 0.03 is in ex- 
cellent quantitative agreement with microscopic model of 
strangeness production and equilibration in QGP [43.l45j. 
adopting latest strange quark mass value. The high QGP 
s-yield oversupplies in hadronization the hadron phase 



15 



space resulting in 7s ~ 2 seen in figure 13 Considering 
the RHIC results shown in figure [151 we see a slightly 
higher s/S saturation limit for most central collisions, 
though the difference is within the RHIC error band (not 
shown). It is possible that s/S LHC2760 resuh is 5-10% 
diluted due to inadvertent inclusion in the entropy count 
of the charm decay hadrons. It is also of interest to note 
that at RHIC62, s/S increases monotonically (discount- 
ing the low confidence level most peripheral point) with 
increasing A'part suggesting that the QGP source reaches 
chemical equilibrium only for most central collisions. At 
LHC2760 there is such increase for much lower size vol- 
ume of the collision centrality A'part < 150. 

In the bottom panel of figure \T5\ we show the ther- 
mal energy cost to make a strange quark-anti-quark pair. 
At LHC2760, the energy cost to make a strange pair is 
practically constant for the wide range of mid-central to 
central collisions, which confirms that strangeness in the 
QGP fireball is in chemical equilibrium at the time of 
hadronization. The slight increase of the thermal energy 
cost for small centralities corresponds to the lower yield of 
strangeness seen in top panel of figure[Tni At RHIC62, we 
see monotonically improving energy efficiency converging 
to a value slightly below our new LHC2760 result, but 
well within the error bar at RHIC62 (not shown). The 
rise of energy cost for smaller systems relates to the fact 
that a larger and notable fraction of strangeness was pro- 
duced in first hard collision processes during the initial 
stages of the collision which for RHIC62 and LHC2760 is 
resulting in higher energy needed to produce one strange- 
anti-strange pair. 



D. Connection to lattice results and related 
considerations 

Elaborate lattice-QCD numerical computations of 
QGP-hadron transition regime are available today [23 . 
[25| . and are comprehensively reviewed in Ref. [2g: the 
HotQCD collaboration f25'| converged for 2-1-1 flavors to- 
wards Tc = 154 ±9 MeV. The question how low the value 
of Tc can be, remains in current intense discussion, as 
the latest work of Wuppertal-Budapest collaboration [27'| 
suggests a low Tc ~ 145 MeV. For an expanding QGP 
with supercooling, this can lead to hadronization below 
Tc ~ 145 MeV and near T = 140 MeV. This is indeed 
the range of values of T that we find in our chemical 
non-equilibrium SHM analysis. 

A comparison of lattice results with freeze-out con- 
ditions is shown in figure 1161 The two bands near to 
the temperature axis display the lattice critical temper- 
ature in the range Tc = 154 ± 9 MeV ^ (red online) 
and Tc = 147 ± 5 MeV [l^l (green online). The symbols 
show the results of hadronization analysis in the T-fis 
plane. We selected here the results for the most central 
collisions and heaviest nuclei. The solid (blue) circles 
are SHARE chemical non-equilibrium results obtained by 
our group, with result presented in this paper included 
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FIG. 16: (color online) Phase diagram showing current lattice 
value of critical temperature Tc calculated by two groups [H, 
[27l |. and results of this work as well as our previous results 
(blue circles) [37i . [4^ and results of other groups [i,[ii,|43- 
152 ]. Full circles refer to chemical non-equilibrium, all other 
symbols refer to fit results with chemical equilibrium of light 
quarks. 



in the LHC domain, and RHIC and SPS results seen, e.g., 
in [nl 113, Eil. The LHC2760 freeze-out temperature is 
in our case clearly below the lattice critical temperature 
Tc- As just discussed, this is expected for supercooling 
followed by sudden hadronization. We show also — 1 
resuhs of other groups: GSI [13, Ell, Florence [H El, 113, 
THERMUS [m, STAR [H and ALICE These results 
show the chemical freeze-out temperature in numerous 
cases well above the lattice critical temperature Tc, which 
in essence means that these SHM calculations are incom- 
patible with lattice calculations. 

The two recent lattice results shown in figure [16] chal- 
lenge the hadronization dogma of the past decade which 
has been adopted without verification due to simplic- 
ity of numerical execution: this ad-hoc GSI-chemical- 
equilibrium-hadronization [l^ produces a hadronization 
temperature above lattice phase cross-over results. Two 
conspiring hypotheses were made in Ref. [jU|: 1) there 
is chemical equilibrium in: 2) a long lived hadron gas 
phase. Both statements were assumptions without the- 
oretical or experimental evidence 'confirmed' by fits to 
data, which had even with the large experimental errors 
a rather large and thus a negligible confidence level. 
Therefore, this model needed additional support. Lattice 
results showing Tc = 173 ±8 MeV were often introduced 
in support of equilibrium-SHM. These can be seen in a 
superficial glance in figure 10 in 53 1. 

Reading the paper [U, one sees that the high Tc — 
173 ± 8 MeV applies to the mathematical case of two 
light quark flavors on discrete space-time. Allowing for 
strangeness flavor in QGP, the hadronization tempera- 
ture must decrease. Therefore, already a decade ago 
Tc — 154 ± 8 MeV was the best estimate for 2-1-1 flavors. 
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leading to the consensus range Tc — 163 ± 15 MeV before 
continuum limit. Present day continuum value we esti- 
mate to be Tc ~ 150 ± 7 MeV combining the two results 
seen in figure fT6l However, many experts w ill p refer the 
more precise lattice with Tc = 147 ± 5 MeV [Ij . To con- 
clude our argument, the hadron gas equilibrium model is 
strongly disfavored in the LHC and RHIC reaction energy 
range by a poor data fit relative to non-equilibrium, and 
by the inconsistency with the lattice results which show 
T > Tc- The chemical non-equilibrium SHM is favored by 
offering an excellent data fit, and the result T < Tc, with 
the difference AT = Tc — T being just where we would 
like it, considering the dynamics of the fireball expansion. 

An important requirement, for the full chemical 
non-equilibrium hadronization approach, is that in the 
hadronization process, quark flavor abundances emerge 
as produced at an earlier and independent stage of fire- 
ball evolution. Our analysis relies on hadronization be- 
ing fast, not allowing a significant modification of the 
available quark abundances. These quark abundances at 
LHC in a wide range of centralities and in most central 
RHIC collisions are near to the QGP chemical equilib- 
rium abundance. In order for the quark yields to remain 
largely unchanged during hadronization and after, it is 
necessary that the transformation from QGP to hadrons 
(hadronization) occurs suddenly and at relatively low 
temperature, near to the expected chemical freeze-out 
point where particle abundances stop evolving. There 
has been an overwhelming experimental evidence favor- 
ing sudden hadronization. This has been seen in pion cor- 
relation (HBT) results for a long time [2l|. The sudden 
hadronization model was required for consistency with 
these results [H, [23| ■ It is associated with chemical non- 
equilibrium SHM analysis of the data [1, @ . Today, with 
lattice QCD transition conditions reaching a low T con- 
sensus the only SHM approach that remains valid is the 
chemical non-equilibrium. 



IV. DISCUSSION 
A. What is new at LHC 

The primary difference between RHIC62 and LHC2760 
data is a 4-times larger transverse volume dV/dy at 
hadronization, as seen in top section of figure |9l This 
result is in qualitative agreement with the pion correla- 
tion HBT studies [H^l • Given the nearly constant entropy 
density at hadronization, the growth of volume drives the 
total entropy yield, which is up to 3.2 times greater at 
LHC2760 than at RHIC62. 

Other differences of LHC2760 compared to RHIC62 are: 
1) an order of magnitude smaller baryochemical poten- 
tial /iB — 1.5 MeV, see figure (TU] 2) constant as function 
of centrality phase-space occupancy 7^; 3) earlier satura- 
tion of 7s as a function of centrality, and thus 7s/7g-ratio 
following the behavior of 7s. For comparison, at RHIC62, 
we have a fast increase of 7s over the entire range of A^part , 



as is shown in bottom panel of figure IHl The LHC2760 
result is interpreted to mean that the QGP fireball is 
rapidly chemically equilibrated already for small iVpart, 
while at RHIC62, we must have a large value of A^part, 
that is a large volume, and thus large lifespan, to achieve 
full strangeness chemical equilibrium in the QGP fireball. 
The value s/S — 0.03 is in excellent qualitative agree- 
ment with microscopic model of strangeness production 
and equilibration in QGP and the associated predictions 
of the final state yield [H, HI] . 

As a comparison of our present work with our pre- 
dictions [3^ shows, the yield of strangeness is ^ 20% 
below our prior expectations. These were motivated by 
consideration of a very rapidly diluting QGP fireball, 
wherein the early strangeness QGP equilibrium is pre- 
served and leads to overabundance, above QGP chem- 
ical equilibrium at time of hadronization. Such a be- 
havior was indicated given the RHIC results showing a 
steady rise, see figure [15] for RHIC62. Instead, we find 
a perfectly equilibrated QGP fireball: the observed value 
of s/S* ~ 0.03 is expected for a chemically equilibrated 
QGP fireball near hadronization condition. This equi- 
librium QGP saturated value s/S = 0.03 is observed for 
many centralities. Since to obtain our prediticon we used 
s/S — 0.037 both the value of 7s and yields of Kaons are 
equally ~ 20% suppressed compared to expectation 's^ , 
and other strange particles as well. How this is possible 
will be one of the riddles that future data and theoretical 
modeling will need to address. For us, this strangeness 
suppression compared to expectation is the most remark- 
able difference from RHIC data that we have found in this 
first LHC result analysis. 



B. Centrality dependence 

Considering the bulk properties of the fireball at 
hadronization, the most remarkable finding is that there 
is so little centrality dependence. This means that at 
LHC2760 the source of hadrons is a hot drop of en- 
ergy that varies mainly in volume as we vary the col- 
lision geometry. This applies to energy density s ~ 
0.50 ± 0.05 GeV/fm'^, hadronization pressure P and en- 
tropy density a in the entire centrality range, see table Hill 
and figure [131 These bulk properties decrease monoton- 
ically and slowly and assume the smallest value for the 
most central collisions, supporting the reaction picture 
of expanding and supercooling fireball — the larger sys- 
tem supercools a bit more. Recall that the error bands 
in figure [T31 are based on 7s uncertainty. The one clear 
centrality dependence of the fireball we find is the rapid 
rise and early appearance of the strangeness yield satu- 
ration seen in figure [T5l 

The chemical freeze-out temperature T decreases by 
about 3 MeV at all centralities compared to RHIC62, see 
middle panel in figure [SI (we do not consider here the 
most peripheral RHIC62 result which has small confi- 
dence level). We believe that this result is related to 
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the need to expand and supercool further the initial 
energy and entropy rich LHC2760 fireball. The large 
expanding QGP matter pushes further out, supercool- 
ing more and yielding a further reduction in the sudden 
hadronization temperature. The freeze-out temperature 
T increases considering more peripheral collisions, see fig- 
ure Eb, which can be explained by the disappearance of 
supercooling present for the most central and most en- 
ergetic collision systems. Considering the behavior of 
both LHC2760 and RHIC62 for TVpart ^ 0, we obtain 
Tliad ^ 145 ± 4MeV, applicable to hadronization with- 
out supercooling. This value is in good agreement with 
the latest lattice result [27[ for transformation tempera- 
ture from QGP to hadrons. 



The value of temperature and its behavior as a func- 
tion of centrality and heavy ion collision energy suggest 
that produced hadrons emerge directly from a sudden 
break-up of quark-gluon plasma. The hadron particle 
density at this low T is sufficiently low to limit the parti- 
cle number changing reactions and render these insignif- 
icant. T = 145-140 MeV is at, and below, the expected 
QGP phase transition. The presence of chemical non- 
equilibrium at this low T means that hadrons did not 
evolve into this condition, but must have been produced 
directly from the deconfined phase. This is consistent 
with the HBT time-parameter, which suggests that par- 
ticles are produced at a scale which is sudden compared 
to the size of the system, as is expected for a super- 
cooled QGP state undergoing, e.g., a filamenting breakup 
at T ~ 140 MeV, and the result of such dynamics is qual- 
itatively consistent with the features described here [21j . 



Last row in table IIIII shows the ratio of entropy at 
LHC2760 to RHIC62, S'lhc/'S'rhic , within the rapidity 
interval —0.5 < y < 0.5. The entropy enhancement fac- 
tor increases monotonically with centrality, from ratio of 
1.27 in most peripheral bin to ratio 3.23 in most central 
bin. This increase requires volume dependent additional 
entropy production mechanisms, which are more effective 
for the more central, larger iVpart, collisions. Such an in- 
crease can arise from hard parton collision generated jets, 
which are better quenched in the larger volume of matter, 
and in addition in abundant charm production, which de- 
cays into hadrons and appears as additional hadron mul- 
tiplicity, i.e., entropy. As long as the additional entropy 
is generated in early stages of the fireball evolution, this 
has little impact on SHM method of approach in study of 
hadronization. For example, the quenching of QCD jets 
feeds thermal degrees of freedom that can convert a part 
of its energy into strangeness. However, charm decay is 
different as it occurs after hadronization. Thus, it needs 
to be accounted for and/or proved irrelevant. It is possi- 
ble that charm decay entropy generating mechanism may 
be the cause of the slight (5%) strangeness s over entropy 
S dilution at LHC2760 (see figured!]). 



C. What we learn about hadronization at LHC 

The full chemical non-equilibrium is introduced by the 
way of the parameter 7^ 7^ 1. This allows to describe 
a situation in which a source of hadrons disintegrates 
faster than the time necessary to re-equilibrate the yield 
of light quarks present. The HBT data provide experi- 
mental evidence that favors a rapid breakup of QGP with 
a short time of hadron production 12111 . and thus favors 
very fast, or sudden, hadronization [2a f23j . There has 
been for more than a decade an animated discussion if 
this parameter is actually needed with arguments such 
as simplicity used to invalidate the full chemical non- 
equilibrium approach. 

We have shown that only the chemical non-equilibrium 
SHM describes very well all available LHC2760 hadron 
production data obtained in a wide range of geometric 
centralities obtained in the rapidity interval —0.5 < < 
0.5, and the outcome is consistent with lattice QCD re- 
sults. We successfully fit the data with x^/ndf < 1 for all 
centrality bins, and show a smooth systematic behavior 
as a function of centrality of both, the statistical SHM 
parameters, see figure [9l and bulk physical properties, 
see figure I13[ that allow a simple and consistent inter- 
pretation. SHM is validated at LHC2760 as it describes 
precisely yields of different particles in a wide range of 
geometric centrality of the collisions and which span over 
more than 5 magnitudes, see figure [51 

We have shown that it is impossible to fit the ra- 
tio p/tt = 0.046 ± 0.003 0, Q together with the other 
data, when choosing a SHM with — 1. However, 
p/tt ~ 0.05 is a natural outcome of our chemical non- 
equilibrium fit where 7, ~ 1.6. This result was pre- 
dicted [3^: within the chemical non-equilibrium SHM, 
p/tt I prediction = 0.047 ± 0.002 for P = 82 ± 5 MeV/fm^ is 
in agreement with experimental result we discussed here, 
for most central collisions p/ttIalice = 0.046 ± 0.003. 
We do not therefore understand the comment by ALICE 
collaboration about this prediction: "actual predictions 
underpredict the measured ratio" 's^. 

We have discussed, in section ITlI Bl the possibility of 
p/tt ratio evolving after hadronization, and found this 
scenario to be highly unlikely considering that experi- 
mental ratio p/tt does not vary in a wide centrality do- 
main. Therefore, the fact that chemical equilibrium SHM 
variant over-predicts p/tt and produces a poor Xtot^iJ see 
figure lll[ demonstrates that the chemical equilibrium 
SHM approach does not work at LHC2760. Further evi- 
dence for the chemical non-equilibrium SHM comes from 
universality of hadronization at LHC2760 and at RHIC, 
see subsection IIII CI and Ref . l38l . 



D. Predicting experimental results 

Our prediction of hadron yields [1^ required as in- 
put the charge particle multiplicity dNch/dy which nor- 
malizes the reaction volume dV/dy. Further, we as- 
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sumed strangeness per entropy content s/S, and the 
nearly universal hadronization pressure with preferred 
value P — 82 zL 5MeV/fm'^. This is accompanied by 
the strangeness conservation constraint (s — s) = and 
the projectile-target charge to baryon ration Q/B = 0.4 
and, as baryocheniical potential cannot yet be fully de- 
fined, an approximate value 0(1) MeV. Using this input 
with a 5% error, we obtain the most compatible values 
of dV/dy, T, 7^, 7^ and chemical potentials, and we can 
evaluate the particle yields along with fireball properties. 

We have redone the predictions for y^s^iy — 2.76 TeV 
case with the tested and released SHAREv2.2 code and 
find that the pre-release SHARE predictions in [s^ were 
made for dNch/dy = 2150 and not for dNch/dy = 1800. 
Therefore, all absolute hadron yields stated in Ref. [s^ 
are normalized to be ~ 20% too large, in addition to 
the strangeness overcount originating in the assumption 
s/S = 0.037 > 0.030. The ratios of hadrons with same 
strangeness content were correctly predicted. 

Applying our prediction method using the updated 
strangeness value oi s/S — 0.030 and a more precise 
hadronization pressure estimate P ~ 77 ± 4, MeV/fm 
results, for ^/snn — 2.76 TcV, in accurate prediction 
of all hadron particle yields, statistical parameters, and 
fireball bulk properties, without using as input any in- 
dividual hadron yield. This validates our approach [39j. 
which can be applied to the forthcoming Pb-Pb colli- 
sions at y^SNN = 5.5 TeV or in the RHIC beam energy 
scan. Noting that the multiplicity of produced hadrons is 
synonymous to entropy of the fireball, this result means 
that all hadron yields can be predicted within the frame- 
work of chemical non-equilibrium SHM using as input the 
properties of the bulk matter in the fireball. 



E. Conclusions and outlook 

The good fit of all observed particles allows us to 
predict with some confidence the yields of yet unmea- 
sured hadrons within the chemical non-equilibrium SHM 
scheme, which are seen in table |lTl The question is how 
stable these yields are when data basis of the fit increases 
to include new measurement. A small SHM parame- 
ter change should be expected also when we refine the 
theoretical model by adding features, such as inclusion 
of hadrons from perturbative QCD jets and/or charm 
hadron decay contribution to hadron yields. We believe 
that predictions for the primary 'stable' hadrons such 
as 77 are accurate. On the other hand, even the minor 
changes in SHM parameters can have relatively large ef- 
fect especially for anti-matter clusters shown at bottom 
of tableini in the anti-alpha, we have 12 anti-quarks, and 
a few % error in understanding their primordial yield is 
raised to 12th power. 

It is quite remarkable that despite a change by a fac- 
tor of 45 in reaction energy, we find for all centralities 
at both LHC2760 and RHIC62, that the energy density 
of hadronizing matter is 0.50 ± 0.05 GeV/im^, as is seen 



in figure [T31 In fact, the present day data favor a sys- 
tematic decrease of hadronization pressure P from pe- 
ripheral towards central collisions as compared to ear- 
lier RHIC62 [13, RHIC200 [l3| and our preliminary LHC 
analysis with limited data set ^3 ■ It is possible that the 
more dynamical expansion of the LHC2760 fireball and 
deeper supercooling of the fireball are the cause. 

We checked that assuming universal hadronization 
pressure, we could obtain a very good fit to particle data 
for all centrality LHC2760 data bins. This means that 
if and when more hadron yield data are available, the 
decrease in bulk properties with centrality seen in fig- 
ure [13] could easily disappear. Therefore, the presence 
of a constant critical hadronization pressure [37| could 
extend from SPS to LHC. We are investigating this hy- 
pothesis, as well as the possibility that another quantity 
governs universality of hadronization. We hope to return 
to the matter as soon as we have understood better the 
final state contributions to hadron yields from charmed 
hadron decays. 

We have shown that the precise hadron yields mea- 
sured by the ALICE collaboration at LHC2760 have of- 
fered a vast new opportunity to explore the properties of 
the QGP fireball and to understand the dynamics of its 
evolution and matter production. We are able to quan- 
tify the key physical properties at this early stage. With 
more data becoming available, we expect a significant re- 
finement and improved understanding of both the QGP 
fireball and mechanisms of matter creation out of the 
deconfined QGP phase. 
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Appendix A: Data rebinning 

Rebinning multi strange hadron yields 

Since there is no literature stating explicitly the yields 
of S and SI in Pb-Pb collisions, we proceed to obtain 
these results by unfolding the preliminary enhancement 
data. We combine the yield of S and fl produced in p-p 
collisions at 7 TeV [34)] stated in table IIVI and labeled 
'pp' in the third column therein, with the 'preliminary' 
enhancement E relative to p~p and normalized to a pair 
of participating nucleons shown in Ref. fsij and which we 
also show in the fifth column of table IIVI We generate 
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TABLE IV: Enhancement of multi-strange baryon yields per participant pair relative to p-p collisions, p-p yields and calculated 
yields in Pb-Pb, which we use as input to our interpolation as a function of A^pait- 



Particle 


Ref. 


Centr. 


(iVpart) 


p-p data [dN/dy)pp and enhancement E 


(dAf/dy)pbPb = £0.8 






L34J 


p-p 


2 


(8.0 ±0.7) X lO-" 










60-90% 


17.6 


1.58 ±0.18 


0.090 ±0.010 






40-60% 


68.8 


£ = 2.48 ±0.26 


0.55 ±0.06 




[32] 


20-40% 


157 


£ = 2.95 ±0.32 


1.51 ±0.17 




0-20% 


308 


£ = 3.08 ±0.33 


3.08 ±0.33 




[34] 


p-p 


2 


(7.8 ±0.7) X 10"^ 










60-90% 


17.6 


E = 1.57 ±0.19 


0.087 ±0.011 






40-60% 


68.8 


£ = 2.56 ±0.26 


0.56 ±0.06 




[32] 


20-40% 


157 


£ = 3.20 ±0.35 


1.59 ±0.17 






0-20% 


308 


£ = 3.00 ± 0.32 


2.91 ±0.32 


n- 


[34] 


p-p 


2 


(0.67 ±0.08) X 10"^ 












(0.68 ±0.08) X 10"^ 


n- 








60-90% 


17.6 


£ = 2.56 ±0.53 


0.012 ±0.003 


0.012 ±0.003 




[32] 


40-60% 


68.8 


£ = 4.57 ±0.79 


0.08 ±0.02 


0.08 ± 0.02 


20-40% 


157 


£ = 5.23 ±0.95 


0.22 ±0.05 


0.21 ±0.04 


2 




0-20% 


308 


£ = 6.97 ± 1.27 


0.57 ±0.12 


0.56 ±0.12 



the first data point for the centrality bin 0-20% by av- 
eraging the number of participants in the centrality bins 
from to 20% shown in table 1 of Ref. [29] . We reduce 
the yields of both S and by a constant factor of 0.8 in 
order to compensate for the difference in collision energy 
-^i = 7TeV in p-p collisions and y^s^vjv — 2.76 TeV in 
Pb-Pb. We obtained the magnitude of this energy cor- 
rection factor by comparing with the actual yield for the 
0-20% centrality bin given in Ref. To disentangle 

the combined yield of i7 + 57, we use the separated il and 
yields from p-p collisions [13] , see table ITVl 
We use the relative errors of the enhancements to esti- 
mate the errors of the multi strange baryon yields, that 
is ~ 11% for S and ~ 20% for 57. Our adopted 57 error 
is larger by ^ 3% than the error of its yield in 0-20% 
centrality bin ^2§\. We adopted this slightly increased 
error to account for procedure which lead us to estimate 
the yield of 57, 57 based in part on 57 + 57 yield. The 
mathematical operations leading to the yields, the yields 
and widths we use are stated in self-exlanatory fashion 
in table lYl 

To account for the different centrality bins for multi 
strange baryons as compared to tt, K, p and 4>/ K, we 
express the centrality bins in terms of average number of 
participants according to p9l | and then interpolate every 
particle yield dN/dy available as a function of A'part with 
a power law 



dN (A^pa,t) 
dy 



c, 



(Al) 



where a,b and c are free parameters. The form of the 
function has no immediate physical motivation, it serves 
well the purpose of unifying the data across incompatible 
centrality bins. This method enables us to evaluate the 
invariant yields for any given A'part, i-C-, arbitrary cen- 
trality and thus enables us to include the multi-strange 
baryon yields in this analysis. Interpolation parameters 



together with of each particle interpolation are sum- 
marized in table [V] For completeness, and potential fu- 
ture use, we present also the parametrization of tt^, 
and p^ which do not require rebinning. Small values of 
X^/ndf show that our description is accurate in the given 
interval of iVpart ■ Interpolation curves are depicted with 
dashed lines in figure IHl for particle yields. 



TABLE V: Interpolation parameters of particle yields as de- 
fined in Eg. ED 



Particle 


a 


b 


c 


xVndf 


^+ 


0.725 


1.160 


-0.890 


0.30/7 


7r~ 


0.724 


1.160 


-0.864 


0.22/7 


K+ 


0.0935 


1.187 


-0.174 


0.07/7 


K" 


0.0927 


1.188 


-0.158 


0.06/7 


P 


0.0432 


1.116 


-0.047 


0.32/7 


P 


0.0502 


1.089 


-0.088 


0.26/7 




0.00552 


1.108 


-0.043 


0.12/1 




0.00762 


1.049 


-0.067 


0.70/1 


Q- 


0.000286 


1.324 


-0.0006 


0.40/1 




0.000309 


1.308 


-0.0011 


0.21/1 


K*7K 
A/tt 


see text for details 
see text for details 


0.032/1 
0.0054/1 



Rebinning K''*/K , A/tt hadron ratios 

We include particle ratios K°*/K~, (/)/K~ and A/tt = 
2A/(7r- +7r+) [3l|. This adds A, K*° and (j) into our data 
set. Since in some ratios certain systematic uncertainties 
of individual yields cancel out, introduction of ratios is 
reducing the overall error of the global fit. The ratio 
0/K has an experimental data point in each centrality 
bin we analyze, removing the need for interpolation or 
rebinning. Thus the following only addresses K''*/K~, 
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and A/tt. 



TABLE VI: Experimentally measured ratios used as input to 
our interpolation as a function of Apart- 



Ratio 


Ref. 


Centr. 


(Apart) 


Experimental ratio 






60-80% 


22.6 


0.333 ± 0.084 


k*7k- 




40-60% 


68.8 


0.285 ± 0.061 


20-40% 


157 


0.245 ± 0.066 






0-20% 


308 


0.194 ±0.051 






60-80% 


22.6 


0.0355 ± 0.0041 


2A 




40-60% 


68.8 


0.0371 ± 0.0042 


mi 


20-40% 
10-20% 


157 
261 


0.0365 ± 0.0042 
0.0355 ± 0.0041 


(7r-+ 7r+) 






0-10% 


357 


0.0336 ± 0.0040 



There are four, resp. five, data points for K*/K , resp. 
A/tt, which we present in table IVTl We describe K*/K 
dependence on iVpart with a power law 

^=f (^patt) = 16.23( V,,°"34 _ 1) + 0.512, (A2) 

with total x^/ndf = 0.032/1. Systematic behavior of 
A/tt as a function of centrality is qualitatively different 
from K*/K, see figure [3 a power law is not sufficient to 
properly describe the data. We use a sum of two power 



laws in the following form 

- =/(A^part) 
TT 

= -6.79 X IQ-^ 7V°^^.f - 2 Afpa°r"^ + 2.03 , (A3) 

with x^/ndf = 0.0054/1. In these two cases, the form of 
the ratio functions has no immediate physical meaning, 
it is invented in order to provide an accurate empirical 
description; note that the bottom of the table IVl presents 
the fit quality of these two ratios. 

Interpolation curves are depicted with dashed lines, in 
figure[7l for K*/K and A/tt ratios. To assign an error to 
the interpolated data points, we take the average nearby 
experimental error for the given particle yield or ratio. 
However, by extrapolating the K*/K ratio to iVpait — 
382, we introduce systematic error due to our choice of 
the functional form of Ea. E2] To account for this effect, 
we multiply the error of K*/K by 2 (resp. 1.5) in the 
most (resp. second to most) central bin we analyze as 
indicated by the shaded area in figure [71 As seen in 
figure ini we also extrapolate but we do not believe 
that this adds to the already significant error, considering 
that our power law interpolation functions describe other 
hadron yields up to A^pait — 382. 
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